Orientation selection in lamellar phases by oscillatory shears 
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In order to address the selection mechanism that is responsible for the unique lamellar orientation 
observed in block copolymers under oscillatory shears, we use a constitutive law for the dissipa- 
tive part of the stress tensor that respects the uniaxial symmetry of a lamellar phase. An interface 
separating two domains oriented parallel and perpendicular to the shear is shown to be hydrodynam- 
ically unstable, a situation analogous to the thin layer instability of stratified fluids under shear. 
The resulting secondary flows break the degeneracy between parallel and perpendicular lamellar 
orientation, leading to a preferred perpendicular orientation in certain ranges of parameters of the 
polymer and of the shear. 
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Oscillatory shears are often used to promote long range 
order in lamellar phases of block copolymers, yet the 
mechanisms responsible for selecting a particular lamel- 
lar orientation relative to the shear remain unknown. 
A possible mechanism based on a hydrodynamic insta- 
bility in a microphase separated copolymer is presented 
here that can distinguish between the experimentally im- 
portant cases of parallel and perpendicular orientations 
(Fig. ^| . The instability occurs at the boundary separat- 
ing parallel and perpendicular regions provided that the 
dissipative part of the stress tensor of the copolymer is 
chosen to reflect the uniaxial symmetry of these broken 
symmetry phases. Our results rely solely on the uniax- 
ial symmetry of the microphases, and therefore should 
generally apply to other complex fluids of the same sym- 
metry. 

Block copolymers are being extensively investigated as 
nanoscale templates for a wide variety of applications 
that include nanolithography 0, 0, 0] i photonic compo- 
nents 4], or high density storage systems However, 
given the small wavelength of the microphases (tens or 
hundreds of Angstroms) , macroscopic size samples do not 
completely order through spontaneous self assembly. In- 
stead, oscillatory shears are commonly introduced in or- 
der to accelerate long range order development over the 
required distances (see Ref. for a recent review). In 
practice, a variety of lamellar orientations are observed 
depending on the architecture of the block and the pa- 
rameters of the shear 0, IS IS EH El , while the mech- 
anisms responsible for orientation selection are not yet 
understood. The first theoretical analysis of orientation 
selection in block copolymers was conducted in the vicin- 
ity of the order-disorder transition of the copolymer, and 
addressed the effect of a steady shear on the growth of 
critical fluctuations 0. Fluctuations along the perpen- 
dicular orientation were shown to be less suppressed by 
the shear, and hence it was argued that this orienta- 
tion would be selected. Consideration was later given to 
anisotropic viscosities of the microphases, which led to 



different relative stabilities of uniform parallel and per- 
pendicular configurations due to the different effects of 
thermal fluctuations on each orientation 01 ■ Later work 
focused on the role played by viscosity contrast between 
the polymer blocks [14j , and showed that the perpendicu- 
lar alignment dominates for high shear rates, and parallel 
otherwise. Existing experimental phenomenology con- 
cerning orientation selection is far more complex than 
these analyses would suggest, and is seen to drastically 
depend on the architecture of the block jfj . The analysis 
that we present does not rely on fluctuation effects near 
critical points, allows for oscillatory shears, and explicitly 
incorporates hydrodynamic effects resulting from viscos- 
ity contrast between the microphases, thus overcoming 
the limitations of previous treatments. 

The experimentally relevant range of shear frequen- 
cies is well below the inverse characteristic relaxation 
times of the polymer chains, and hence a reduced descrip- 
tion in terms of the monomer volume fraction is adopted 
[lH Hrl flfij |. According to this description, the lamel- 
lar phase response is solid like or elastic for perturba- 
tions directed along the lamellae normal, and fluid like 
or viscous on the lamellar plane. In the limit of vanish- 
ing frequency, the viscous part of the response has been 
assumed to be Newtonian with uniform shear viscosity, 
and therefore parallel and perpendicular orientations are 
degenerate and unmodified by the shear. We address be- 
low the consequences of what we believe is the leading 
deviation away from Newtonian response in the limits of 
low frequency and characteristic flow scale much longer 
than the lamellar spacing: the viscous stress tensor of a 
lamellar phase is, by reason of symmetry, the same as 
that of any other uniaxial phase (e.g., a nematic liquid 
crystal [171 Il8|). The slowly varying local wavevector of 
the lamellae plays a role analogous to that of the direc- 
tor in a nematic. The rest of this paper is devoted to 
the study of the effect of this assumption on the hydro- 
dynamic stability of the configuration shown in Fig. ^ 

We assume that the dissipative part of the linear stress 
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FIG. 1: Schematic representation of a parallel/perpendicular 
configuration begin sheared along y. 



tensor is that of a uniaxial, incompressible phase (18| . 

a B = r/D, tj + aihihjhkhiDki + a 56 (nihkD jk + hjh k D ik ) 

(1) 

with i, j, k — x,y, z, Dij = diVj + djVi with Vi the local 
velocity field, and n = (n x ,h y ,n z ) denoting the slowly 
varying normal to the lamellar planes. The Newtonian 
viscosity is 77, and ot\ and a§§ are two independent vis- 
cosity coefficients. The dynamic viscosity 77' (77' = G" /u, 
with G" the loss modulus and lo the shear frequency) 
is 7/ = 77 for a fully ordered perpendicular configura- 
tion (n = (1,0,0)), and 77' = 77 + for a parallel ori- 
entation (n = (0,0,1)). This model is consistent with 
low frequency rheology in PEP-PEE diblocks showing 
^par > ^perpJuJ if a 56 > 0, and with PS-PI diblocks 
Vpur < Vpeip UM if a 56 < 0- Given this assumption, the 
effective dynamic viscosity in the two domain configu- 
ration of Fig. ^ would be different in each domain. It 
is known that the analogous configuration for the case 
of two Newtonian fluids of different viscosity is unstable 
both for steady 0, H3| and oscillatory [2l| shears. We 
describe below the extension of these results to uniaxial 
phases in the limit of small but nonzero Reynolds number 
flows. 

We consider a base state involving perpendicular (A) 
and parallel (B) regions separated by a planar surface, 
subjected to an imposed shear Vo = jlooI cos(ujt)y at 
z = d, and Vo = at z — (Fig. QJ, with 7 the 
shear amplitude and d the system thickness. The re- 
sulting velocity field va,b = (0, Va,b,0) is the same as 
that of two superposed Newtonian fluids with viscosi- 
ties /j,a = f\ and /is = 77 + CV56. We then consider a 
small perturbation of the domain boundary as shown 
in Fig. ^ and write the velocity fields in A and B as 
vf' B = V A ,B&iy + uf ' B (i = x,y,z). By expanding 
uf' B = £ gx ,g uf' B (qx,q y ,z,t)ex-p[i(q x x + q v y)}, sub- 
stituting into the modified Navier Stokes equation that 
results from the choice of Eq. (I}, linearizing, and elimi- 



nating pressure and u y , we find 

Re [(d t + iq y V A ){dl - q 2 )u* - iq y (d 2 z V A )ui] 
= {dl - q 2 ) 2 uf - a 56 ql(d 2 z - q 2 )uf 
+iq x [2a iq 2 x - a 56 (d 2 ~ q 2 )]d z u^, (2) 
Re [d t (d 2 - q 2 )u^ + iq y (d 2 z - q 2 )(V A U^) + 
2q x q v {d z V A )uf} = (1 + a 56 )(8 2 - q 2 ) 2 u* 

~2a iq 2 (d 2 - q 2 y )u^ (3) 

for the perpendicular domain A (0 < z < d^). Here 
q 2 = Qx + q 2 arm ^ e = pud 2 /r] is the Reynolds number, 
with p the copolymer density. Similarly 

Re [{d t + iq y V B )(d 2 z - q 2 )u B - iq y (d 2 V B )uf] 

= (1 + a 56 )(d 2 - q 2 ) 2 u B - 2a ig 2 d 2 u B , (4) 
Re [8 t (d 2 - q 2 )u B + iq y (d* - q 2 )(V B u B ) + 
2q x q y (d z V B )uf] = (d 2 z - q 2 ) 2 u B 
+a 56 (d 2 - q 2 ) (d 2 u B - iq x d z u B ) - 2 l q x a 1 d 3 z uf , (5) 

for the parallel region B {cLa < z < d). All quantities 
have been made dimensionless by a length scale d, a time 
scale and rescaling viscosities by 77, i.e., a± — > ai/77 
and a 56 — > a 56 /?7 (ha = 1, Pb = 1 + atm)- Rigid bound- 
ary conditions are used on the planes z = and z = d 
(d = 1 after rescaling). At the interface z — dA+C,{x, y, t) 
we have v A = v B , and (a^-o$)fy = — r'(Sg +flg)C*W*i* 
with r' = T/(i]ujd) and T the interfacial tension. Also, 
the kinematic boundary condition for the interface is 

(a t + v s • v)c = vf- 

Equations are similar to the Orr-Sommerfeld 

equation for the Newtonian case 0, I20I Elj . except 
that x and z velocity fields are now coupled (except 
at q x = 0). The solution can be found by writing 
uf: B — exp((T7j)0^^f and £ = exp(at)h, with £ the 
Fourier transform of C, and a the Floquet exponent yield- 
ing the perturbation growth rate (coefficients in Eqs. 
|JSJ proportional to Va,s are periodic in time). 

For typical block copolymers, p ~ 1 g cm~ 3 , ~ 1 cm, 
and 77 ~ 10 4 -10 6 P, resulting in Re/u = 10" 4 -10" 6 s. 
Hence Re <C 1 for the frequencies of interest. We further 
expand the velocity, interfacial functions 4> A,B ' , h, as well 
as the Floquet exponent a in powers of Re, and solve Eqs. 
@-© order by order. In the limit Re — > while keeping 
the surface tension V finite, we find o~ = f z0 (q x ,q y )T' , 
with < for all wave numbers q x and q y . Hence the 
interface is stable, indicating coexistence of parallel and 
perpendicular orientations in this limit. 

The situation is different for small but finite values of 
Re. In this case we need to address the order of V as 
well. Typical values of T ~ 1 dyne/cm lead to T'oj = 
10- 4 -10~ 6 s _1 . Given that Re/u = 10" 4 ~10" 6 s, and 
that to ~ Is -1 in typical experiments, we consider the 
distinguished limit T' = 0(Re). Writing V = TxRe, and 
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FIG. 2: Growth rate a/Re as a function of wave numbers q x 
and q y , for ai = 1, «56 = —0.9, 7 = 1, Re = 5 x 10~ 4 , and 
g!a = 1/2. The maximum growth rate is found at q™ ax = 0, 
q™ ax = ±0.61. 
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FIG. 3: Orientation selection for uniaxial systems under os- 
cillatory shear as a function of the layer thickness ratio and 
viscosity contrast. (U, a > 0) denotes the range in which the 
parallel orientation is unstable, leading to the selection of the 
perpendicular orientation. 



(7 = <TiRe, we find 

cri = /Jute, 8»)ri + ~<*V/fi(<7*, gy). (6) 

Here <5 = toAo is proportional to the viscosity contrast 
to = ha/ Lib, with Ao = (cLa + mds)^ 1 , d-B = d— <1a, and 
the function ffi(q x ,qy) depends on the system parame- 
ters cvi, a56, and on d,A, but not on the shear parameters 
7 and w. Whereas ff is always negative, can be pos- 
itive so that Eq. © illustrates the competition between 
the stabilizing effect of surface tension and the destabiliz- 
ing effect of the imposed shear flow. Note that T\ can be 
written as Y\ = 1/We = 9uj~ 2 , with We the Weber num- 
ber and 8 — T/(pd 3 ). Thus, we have a\ — <J\{j 2 ,o->~ 2 ), 
and instability is increased at large shear amplitudes and 
frequencies. 

Typical results for a as a function of wave vector are 
shown in Fig. [2J Unstable wave vectors are near q x = 0. 
This indicates the absence of interfacial modulation along 
the x direction (the transverse direction for perpendicular 
lamellae), and u x — for the associated velocity pertur- 
bations. We have repeated the calculations for different 
ranges of parameters, and found similar results for both 
a and velocity perturbations. 

Based on these results, further progress in determining 
the stability boundary can be made by examining only 
long waves along the y direction (with q = q y ). Equation 
© can be expanded as f^ = —foq 4 + 0(q 6 ), ff x = 
fiq 2 + !iq A + @{q 6 ): so that the Floquet exponent can be 
rewritten as 

<ri = \5 2 ! 2 hq 2 (OfooJ- 2 \^ 2 h)q\ (7) 

where fa > always, as noted above. Both /1 = 
/i(ck56, dA) and /a = f2(&i, &56, <1a) are complicated but 
known functions of their arguments. For small q stability 
is determined by the sign of f± which depends on a^e and 



layer thickness dA, but is independent of shear parame- 
ters. The calculated stability diagram in the {dA/ds, A's) 
plane (/its = 1 + a^) is shown in Fig. [3] Note the sym- 
metry under dA/ds — > (cU/ds) -1 and \xb — > H~b which 
suggests that this instability is related to the known two 
fluid instability produced by viscosity stratification |20| : 
instability occurs when the thinner domain is more vis- 
cous. In the present case, however, viscosity contrast be- 
tween the domains is not caused by fluid stratification, 
but rather because of the effective viscosity contrast be- 
tween lamellae of different orientations. 

The stability analysis discussed thus far is purely of hy- 
drodynamic nature, but can we used to argue that growth 
of unstable modes above threshold leads to orientation 
selection. As can be seen from Fig. ^ parallel lamellae 
are marginal to velocity fields along x and y directions, 
but will be compressed or expanded by nonuniform flows 
along the z-direction. Conversely, perpendicular lamel- 
lae will be distorted by nonuniform flows along x, but not 
by flows along either y or z. Since the instability mode 
is dominated by a velocity field along z, it will lead to 
weak and oscillatory compression and expansion of the 
parallel region, while leaving the perpendicular lamel- 
lae unaffected. The response of an interface separating 
two lamellar phases, subjected to periodic expansion and 
compression and the other marginal, has already been 
addressed in Ref. H2- We showed that the overall free 
energy of the system is reduced by the motion of the 
interface towards the distorted phase (which is storing 
elastic energy during each cycle of the shear) or, in the 
present case, toward the parallel region. In summary, 
Fig. |2| shows the regions of parameters in which paral- 
lel and perpendicular layers would coexist, and those in 
which the layer of perpendicular orientation would grow 
at the expense of the parallel layer. It would be interest- 
ing to test our predictions by examining a system com- 
prising only two layers of different orientations and with 
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varying ratios cLa/Ae to directly address the stability of 
the configuration under shear, and to indirectly measure 
the coefficient from the location of the instability 
threshold. 

Experiments addressing orientation selection always 
involve coarsening of polycrystalline samples with a dis- 
tribution of grain sizes, and hence a range of ratios 
d-A/ds, and a distribution of orientations. It is generally 
argued that lamellar domains with local wavevector not 
on the parallel/perpendicular plane will be eliminated 
from the distribution rather quickly, and hence that the 
selection of a final orientation will be determined by the 
competition between parallel and perpendicular domains. 
Generally speaking, our results imply selection of the 
perpendicular orientation for finite shear frequencies and 
«56 > 0, the latter case appropriate for PEP-PEE but 
not PS-PI blocks if our assumption in Eq. (JTJ holds. If 
«56 < 0, Fig. 0would suggest that a smaller than average 
layer of perpendicular orientation first grows at the ex- 
pense of neighboring parallel layers. Following this initial 
coarsening in which rf^/rfs increases in time, the stability 
boundary in this figure would be reached. It is difficult 
to assess in this case the impact of other dynamical fac- 
tors affecting coarsening such as the effect of an already 
moving boundary and the concomitant flows, or spatial 
correlations built into the distribution of orientations fol- 
lowing this intermediate coarsening [23|. In summary, if 
Q56 < 0, as well as in the limit of Re — > 0, our study indi- 
cates coexistence of parallel and perpendicular domains, 
not inconsistently with experiments addressing orienta- 
tion selection that indicate dependence on processing his- 
tory or experimental details |a iSs 9] (e.g., quenched or 
annealed history of the sample, and the starting time of 
shear alignment). 

Once within the region of instability, it is instructive 
to analyze the dependence of the growth rate of the most 
unstable perturbation on the parameters of the shear. 
This growth rate is given by cr™ ax = (<5 4 f 2 /16 f')7 4 cj 2 , 
with /' = Ofo — S 2 f2j 2 oj 2 /2, and the corresponding most 
unstable wavenumber is g max = (5(/i//') 1 / 2 7w/2, both of 
which increase with shear amplitude and frequency. By 
noting that cr™ ax = a max /Re — a max r]/(pd 2 ui), and that 
usually cr max <C (pd 2 S 2 f 2 /8\f2\r])j 2 uJ for small enough 
0max, we find the maximum growth rate to be given by 
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Therefore, perturbation growth for a given block copoly- 
mer is constant along the line 7W 3 / 4 =const. To our 
knowledge, the only experimental determination of the 
boundary in parameter space separating regions in which 
parallel or perpendicular lamellae are selected has been 
given for PS-PI copolymers 0, [HJ ■ From a limited data 
set, it was approximated by jlj — const. Since it is not 
inconceivable that the experimentally determined bound- 
ary does not correspond to the true stability boundary, 



but rather to the line in which <7 ma x becomes experimen- 
tally observable |24J , it would be desirable to conduct the 
experiment in a block copolymer with > 0. 

In summary, by assuming that the dissipative part of 
the linear stress tensor of a block copolymer has to re- 
spect the broken symmetry of uniaxial lamellar phases, 
we have obtained a long wavelength hydrodynamic insta- 
bility of the interface separating lamellae of parallel and 
perpendicular orientations under an imposed oscillatory 
shear. The instability leads to nonuniform secondary 
flows, which would favor the perpendicular orientation 
in large regions of parameter space. Since our results 
follow from the symmetry of the microphases, we would 
expect them to hold in other complex fluids of the same 
symmetry. 
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